Flavor-Mixing Effects on the QCD Phase Diagram at non-vanishing 
Isospin Chemical Potential: One or Two Phase Transitions? 
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Abstract 



We investigate effects of a fixed nonzero isospin chemical potential on the 
s * ■ T phase diagram of strongly interacting matter using a Nambu-Jona-Lasinio-type 

four fermion interaction. We focus on the influence of a flavor-mixing interaction 
(f-) ■ induced by instantons. We find that already for rather moderate values of the 

coupling strength in the flavor-mixing channel the recent findings of two seperate 
£^ \ phase transitions do not persist. 



The structure of the phase diagram of QCD at nonzero temperature and baryon chem- 
ical potential has been intensively studied throughout the last decade, as well on theoret- 



<]} [ ical as on experimental side. In the baryon chemical potential-temperature (/i#-T) plane 

roughly speaking three domains can be distuingished: the hadronic phase at low tem- 
perature and density, a quark matter phase at low temperature and high baryon density, 
^ which is most probably a color superconductor JTJ |2] , and the quark gluon plasma at high 

temperature. For two flavors at zero chemical potential the hadronic phase is expected 
to be seperated by a crossover from the quark-gluon-plasma 3\. The fact that at zero 
temperature a first-order transition to quark matter is expected [4 implies the existence 
of a second-order endpoint somewhere in the /le-T-plane. This second-order endpoint 
is one of the features of the QCD phase diagram which could be detected in heavy-ion 
collisions jH]. 

Up to now most theoretical investigations of the phase diagram and are restricted to 
zero isospin chemical potential /// (see, e.g., Refs.jTJEJIHE])- The influence of a non-zero 
/ij has been studied in more detail in the context of quark matter at high baryon density 
and low temperature, e.g., [QUI IE], as well as for vanishing \xb P ITU]. 

Recently the [Ab-T phase diagram for a fixed fii ^ has been studied using a random 
matrix model ^JJ and a Nambu-Jona-Lasinio (NJL) model ^2j- The authors of Refs. fTJ 
IT2] find striking effects: There are two first-order phase transitions at low temperature 
and high baryon chemical potential and thus two second-order endpoints. This could be 
important, e.g., for heavy-ion collisions where fij is supposed to be nonzero. In this letter 
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we will argue that the existence of two separate phase transitions becomes unlikely, once 
flavor- mixing effects due to instantons are taken into account. 
Our starting point is the following Lagrangian for two flavors: 

£ = £ + £ 1 + £ 2 , (1) 

with a free part 

£ = q(i$-m)q, (2) 
and two different interaction parts (see, e.g., (HUE]), 

Cl = G^qqf + {qfqf + {qi^qf + {qi^rq) 2 } (3) 

and 

£ 2 = G 2 [{qq) 2 - {qrqf - {qi lb qf + (qi^q) 2 } , (4) 

where q = (u, d) T . G\ and G 2 are coupling contants of dimension energy -2 and m = 
diag{m u ,m l j) contains the current quark masses. For simplicity we assume m„ = m^. 
The interaction is invariant under SUl(2) x SUr(2) x £/y(l) transformations. £i exhibits 
an additional Ua(X) symmetry, whereas this is not true for £ 2 . £ 2 has the structure of a 
't-Hooft determinant in flavor space [T^]. This interaction can be interpreted as induced 
by instantons and reflects the [/^(lO-anomaly of QCD. The reason for this particular 
choice of the two interaction parts will become clear below. 

We now want to study the properties of the system at temperature T and the up and 
down quark chemical potentials 

= yU + <fy« > = ^ - Sfi , (5) 

which are in general different. Here \i = is the quark number chemical potential 

and 5 fx = fii/2 is proportional to the isospin chemical potential *. To obtain the mean 
field thermodynamic potential Q(T, [i, dfi), we linearize the Lagrangian in the presence of 
the following quark condensates: 

(j) u = (uu) , (j) d = (dd) , (6) 

which a priori can be different. We do not consider p = ^((ujsd) — (d^u)), i.e., we 
exclude the possibility of pion condensation. This restricts our model to values of \8fi\ < 
m 7r /2 ~ 70 MeV. For the moment we also do not consider color superconducting phases. 
This will, however, not change our results qualitatively. 

For our further proceeding it is convenient to introduce constituent quark masses: 

M l = m i -4G 1 (j) i -4G 2 (f) j , i ^ j e {u, d} . (7) 

*This follows from the definitions of the baryon density, ub = \{n u + rid), and the isospin density 
n i = \{riu — rid) together with the thermodynamic relation = —dfl/d/u,i for i = u,d,B,I. Note that 
often different definitions are used. For instance, \xb and \ii of Refs. |5lll2| correspond to our fi and 8ti. 
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Using this definition and Eq. (jSJ) for the chemical potentials of up and down quarks, 
respectively, we obtain for the mean field thermodynamic potential: 

fi(T,^,/i d )= ^O (T,/i / ;M / ) + 2G' 1 (^ + ^)+4G 2 ^ (i , (8) 

f=u,d 

where ft (T, fif, Mf) corresponds to the contribution of a gas of quasiparticles of flavor /: 
Q (T, fi f ; M f ) = ~^J dpp 2 [E f + T In (l + exp (~(E f - ///))) 

+ T In (l+exp(~(E f + ///)))} (9) 

with Ef = ^Mj + p 2 . In order to determine the physical solutions, we have to look for the 

stationary points of the thermodynamic potential with respect to the two condensates (j) u 
and 4>d- This leads to the following gap equations which have to be solved selfonsistently: 

fa = ~ J dPP'^ji 1 - < E f) ~ - ( 1Q ) 

where n(Ef) = l/(exp((Ef — jif)/T) + 1) and n = l/(exp((Ef + /i/)/T) + 1) are Fermi 
occupation numbers. Here we have to keep in mind that via Eq. (|7J), the constituent mass 
Mj for one flavor depends in general on both condensates and therefore the two flavors are 
coupled. At this point our separation of the interaction part becomes clear: If we for the 
moment neglect £ 2 , i-e., G 2 = 0, the two flavors decouple. That means Mj only depends 
on the condensate of the same flavor 0j and the mixed contribution to Q (last term of 
Eq. ©) vanishes. In this limit we recover the expression for the thermodynamic potential 
of Ref. ^2] fo r P — 0. In the opposite limit, i.e., G\ = 0, we have "maximal" mixing: The 
constituent mass of flavor i only depends on the condensate <pj with i ^ j. For G\ = G 2 
we recover the original Lagrangian proposed by Nambu and Jona-Lasinio |16j . implying 
M u = M d . 

To study the effects of flavor mixing, let us now write 

= (1 - a) 6?o , G 2 = aG , (11) 

and calculate the phase diagram for fixed Go but different values of a. The amount of 
flavor mixing is thereby controlled by the particular value of a while the values of the 
vacuum constituent quark masses M vac are kept constant. 

For our numerical studies we use the following set of parameters: m u = = 6 MeV, 
a three-dimensional sharp cutoff A = 590 MeV and Go A 2 = 2.435, corresponding to 
vacuum constituent masses M vac = 400 MeV. With these parameters we obtain reasonable 
values for the pion mass, decay constant and the quark condensate in the vacuum: m, = 
140.2 MeV, U = 92.6 MeV and (uu) = (-241.5 MeV) 3 . 

We begin our discussion with the results at T = 0. In Fig.[T]we display the values of M u 
and Md as functions of the quark number chemical potential \x for fixed 5ji = 30 MeV \ 

^Following common practice (e.g., we take a positive value of 8[i, although for the description 

of heavy-ion collisions Sfj, < would be more appropriate. However, since changing the sign of S/i does 
only interchange the roles of up and down quarks, this does not alter our conclusions. 
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Figure 1: Constituent quark masses M u (solid) and Md (dashed) at T = as functions of 
quark number chemical potential \i for 5/j, = 30 MeV and a = (left), a = 0.05 (center), 
and a = 0.11 (right). 



The left panel corresponds to a = 0, i.e., to the case without flavor mixing. We observe 
two distinct phase transitions at fi = 353 MeV for the up quarks and at fi = 413 MeV 
for the down quarks. This behavior is easily understood when we recall that at a = the 
up and down quark contributions to the thermodynamic potential completely decouple. 
Hence, if we had plotted M u and Md, in terms of the corresponding flavor chemical 
potential \i u and fid, respectively, we would have found two identical functions with a 
phase transition at /i/ = 383 MeV. This is basically the result reported in Refs. [TT| I12j. 

Now we want to study the influence of a non-vanishing flavor mixing. In the central 
panel of Fig. show the behavior of the constituent quark masses for a = 0.05. The 
situation remains qualitatively unchanged, i.e., we still find two distinct phase transitions. 
However, because Md now also depends on <ft u (and thus on M u ), and vice versa, both 
constituent masses drop at both critical chemical potentials. Moreover, this small amount 
of flavor mixing already diminshes the difference between the two critical quark number 
chemical potentials considerably. Finally, for a larger than a critical value of 0.104 we 
find only one single first-order phase transition. This is illustrated in the right panel of 
Fig. [H which corresponds to a = 0.11. 

Next, we extend our analysis to non- vanishing temperature. The phase diagrams in 
the fi-T plane for fixed S/i = 30 MeV and three different values of a are shown in Fig. El 
At a = (left panel) we qualitatively reproduce the results discussed in Refs. fH\ IT2]. 
i.e., two separate first-order phase boundaries which end in two second-order endpoints. 
Again, since for a = the up and down quarks decouple completely, we would obtain two 
identical phase diagrams if we plotted the phase structure of flavor / in the fif-T plane. 
In the central panel of Fig. |2]we consider a = 0.11, i.e., slightly larger than the critical 
value a c (T = 0) = 0.104 for a single phase transition at T = 0. Accordingly, there is 
only one phase boundary at low temperatures, but at T = 25 MeV it splits into two lines 
which end at two different second-order endpoints. The two branches are, however, very 
close to each other and already at a = 0.12 we find only one phase boundary with a single 
endpoint. This is illustrated by the diagram on the right, which corresponds to a = 0.15. 
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Figure 2: Phase diagrams in the /i-T-plane for 5/i = 30 MeV and a = (left), a = 0.11 
(center), and a = 0.15 (right). The lines correspond to first-order phase boundaries which 
end in second-order endpoints. 



In our example a rather small amount of flavor mixing is sufficient to remove the 
existence of the second phase transition: Of course, there must be a single phase transition 
at a = 0.5, where M u and M d are equal (see Eq. ((Zj)). (This was the case studied in 
Ref. jH].) However, the critical value a c ~ 0.12 we found in our example is much smaller. 
At T = 0, a rough, but perhaps more general estimate for the critical a can be obtained 
from the observation that the phase transition takes place when the chemical potential 
of quark / comes close to its constituent mass, i.e., /// ~ Mj. Applying this condition 
to the u quark we expect the first phase transition to take place at \i u w M vac , i.e, at 
ji M vac — 5ji. At this point M u drops and, according to Eq. (JJJ) M d drops as well. 
Neglecting the current quark mass, we find 

M d ~ -(1 - a) AG <p d < (1 - a) M vac . (12) 

If this value becomes smaller than the value of fid, we expect also the down quarks to 
exhibit a phase transition. Hence, we estimate 

« C (T = 0) < . (13) 

Note that this estimate would not be affected by a possible restoration of the ^7a(1) 
symmetry at the phase boundary. Obviously, if G2 goes to zero, M d would drop as well. 

For our example, Eq. (|13|) gives a c (T = 0) < 0.15. Comparing this value with the 
numerical result a c (T = 0) = 0.104, we see that Eq. (|13j) is a quite conservative estimate. 
This is easily understood, since in the second step of Eq. (JT2|l we have neglected the fact, 
that (pa also becomes smaller. Our estimate does also not include the observation, that 
the critical chemical potential for the first phase transition rises with a. In any case, we 
have to admit that our arguments cannot explain quantitatively why Eq. (|13|) seems to 
hold even for temperatures approaching the critical endpoint where the quark masses do 
no longer drop discontinuously. 
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At this point one can ask, which value of a is "realistic". To answer this question it 
is helpful to have a look at the 3-flavor NJL model, where the strength of the 't Hooft 
interaction has been fitted by several authors to describe the ^-//-splitting ^3QH|- For 
three flavors, the 't Hooft determinant is a six-point interaction ^H], and the constituent 
quark masses in a 3-flavor NJL model are given by 

M i = m i -4G(f)i + 2K(f) j (f) k , i ^ j ^ k ^ i e {u,d, s} . (14) 

When we compare this with Eq. (JZj) we can identify G\ = G and G 2 = —\K(f) s and thus 

G - ±K(p s 

If we take, for instance, the values of Ref. [Ej, A = 602.3 MeV, GA 2 = 1.835, KA 5 = 
12.36, and <p s = (—257.7 MeV) 3 , we find a ~ 0.21. For the parameters of Ref. [THj we 
get a somewhat smaller value, a ~ 0.16. On the other hand, the success of the instanton 
liquid model to describe hadronic correlation functions [TH] would suggest that C 2 is the 
dominant part of our Lagrangian, i.e., a ~ 1. Anyway, in all cases we would find only one 
phase transition for 5fi = 30 MeV. Typical values of \S/j,\ in heavy-ion collisions are likely 
to be smaller than that. (A simple estimate, assuming the density ratio n u : rid = 290 : 334 
as in 208 Pb, and the approximate relation n u : rid ~ (a*« : /^d) 3 yields 5fi ~ — 10 MeV for 
/j, = 400 MeV. Empirically, one finds fij = Ibii = —5 MeV at chemical freeze-out for 
Pb-Pb collisions at SPS [201 and /x 7 = -12 MeV for Si+Au collisions at AGS jH10 

However, before drawing quantitative conclusions, we should be aware of several short- 
comings of the present model. First, the description of the "hadronic phase" as a gas of 
quarks, rather than hadrons, is certainly unrealistic. In fact, it is quite obvious that 
any prediction of the critical endpoint(s) in non-confining mean-field models should not 
be trusted. (Lattice calculations, although not yet generally accepted, indicate that 
the critical endpoint of QCD at 5fi = is located at \mb = 3// ~ 725 ± 35 MeV and 
T = 160±3.5 MeV [22], i.e., at lower chemical potential and higher temperature than our 
or other NJL results, e.g., ^21E]-) Moreover, as already mentioned, we have neglected 
the possibilities of pion and diquark condensation. Since pion condensation, although 
interesting by itself, does only occur for 5fi > 70 MeV, it is irrelevant for the present 
discussion. A possible diquark condensation in the 2SC phase |T] would further favor a 
single phase transition, since this phase requires an approximately equal number of up and 
down quarks. An extension of the present model to include these possibilities is straight 
forward. 

Keeping the limitations of our model in mind, our results show that flavor-mixing 
effects cannot be neglected in the discussion of the phase diagram. The pure existence 
of these effects is related to instantons and the ?7 J 4(l)-anomaly of QCD. Of course their 
magnitude is a matter of debate, but they probably cancel the interesting phenomena 
discussed in Refs. [TTJ [TJ] . 
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